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Abstract 

The atmospheres of (exo) planets and moons, as well as reflection nebulae, 
contain in general independently scattering particles in random orientation and are 
often supposed to be plane-parallel. Relations are presented for the (bidirectional) 
reflection function and several related functions of such a medium in case the di- 
rections of incidence and reflection both tend to horizontal directions. The results 
are quite general. The medium may be semi-infinite or finite, with or without a re- 
flecting surface underneath, and vertically homogeneous or inhomogeneous. Some 
approximative formulae for the reflection function of a plane-parallel medium with 
independently scattering particles in random orientation, including Lambert's law, 
may be very inaccurate if the directions of incidence and reflection are both nearly 
horizontal. 

keywords: planets and satellites - scattering - radiative transfer - reflection nebulae 

1 Introduction 

A well-known subject in astrophysics concerns multiple scattering of (electro- 
magnetic) radiation in an extended medium containing small, independently scat- 
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tering particles (see e.g. Ambar zumian 1 9431; Chandrasekh ar 1950: Sobol evlll97 



de Hulstll98(]| : lHansen & Travis,1974 : lHovemer et alj2004; Mishchenko et al 



20061) . Examples of such media are provided by the atmospheres of (exo)planets 



and satellites, as well as reflection nebulae and protoplanetary disks. The medium 
is often supposed to be locally plane-parallel, so that one can use the theory devel- 
oped for a plane-parallel atmosphere, i.e. a horizontally homogeneous atmosphere 
of infinite horizontal extent. The radiation (which we will also call light) coming 
from a distant source, like the Sun or a star, may illuminate the top of the atmo- 
sphere and then be scattered by the particles inside before leaving the atmosphere 
at the top in all upward directions. This is called reflected radiation and, neglecting 
polarization, the angular distribution of its (specific) intensity can be expressed by 
means of the so-called reflection function. 

The reflection function is an important fundamental property of an atmosphere, 
normalized so that it is identically equal to one for a perfectly white surface fol- 
lowing Lambert's law. Once the reflection function of an atmosphere has been 
obtained, one readily finds the angular distribution of the intensity of the reflected 
radiation for any angular distribution of incident radiation at the top. In the liter- 
ature this function has a variety of names, such as reflec tion coefficient jSobolevI 
ll975l : lYanovitskiil 19971). bidirection al reflection function dMishchenko et alj 19991) 
and reflectance factor THapkdl 19931) . 

Numerous theoretical and numerical studies of reflection functions have been 
reported in books and papers. Yet, very little attention was given to the limiting 
case when the directions of incident and reflected radiation both tend to horizontal 
directions. This case not only provides more insight into the angular distribu- 
tion of the reflected radiation, but it is also important for exact and approximate 
computations, in particular when discontinuities are involved. This was shown 
for the in tensity of the reflected radiation, first when polarization is negle cted 
faovenier & Stam 2006) and later when it is taken into account (Hovenier & Stanj 
l2007h . As fa r as the refle c tion f unction is concerned, an interesting statement 
was made bv Ivan de Huls namely that the reflection function becomes 

infinitely large for horizontal directions of both incidence and reflection. But he 
restricted himself to the azimuth-independent terms in a Fourier series expansion 
and he did not give any evidence or clarification regarding this statement. 

The principal aim of this paper is to present a comprehensive treatment of 
the behavior of the reflection function and related functions when the directions of 
incidence and reflection both tend to horizontal directions. The organization of this 
paper is as follows. Some basic concepts and definitions are discussed in Sect. [2] 
In Sect. [3] the limiting process of directions of incidence and reflection tending to 
horizontal directions, while keeping the azimuth difference of the directions fixed, 
is considered for, respectively, the reflected intensity, the reflection function and 
orders of scattering of both. Simple examples are given in Sect. |4] for vertically 
homogeneous as well as inhomogeneous media. Section[5]is devoted to functions 
that are related to the reflection function. The azimuth dependence of the incident 
and reflected light is treated in Sect. [6] Approximations are discussed in Sect. |7] 
and some concluding remarks are presented in Sect. [8] 
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2 Concepts and defintions 



We consider a plane-parallel atmosphere composed of randomly oriented particles, 
which may include gas molecules. The particles scatter radiation independently 
and without change of wavelength in all directions with a scattering angle distribu- 
tion called the phase function. Since we are mainly interested in natural particles 
and physically realistic model particles we will assume that the albedo of single 
scattering is positive but not larger than one and the phase function is a positive, 
bounded and continuous function of directions. It is normalized so that its average 
over all directions equals unity. 

There are no internal sources in the atmosphere. A parallel beam of radiation, 
coming from a distant source, is incident on each point of the top of the atmo- 
sphere. The net flux per unit area normal to this beam is ttFq. The direction of the 
incident beam is given by /jq, which is the cosine of the angle this direction makes 
with the downward normal, and an azimuthal angle 0o- Polarization is ignored and 
we focus on the reflected radiation, i.e. the radiation that emerges at the top of the 
atmosphere. Its direction is described by /u, the cosine of the angle this direction 
makes with the upward normal, and an azimuthal angle 0. The azimuthal angles 
are measured from an arbitrary zero direction in an arbitrary sense and only the 
difference, - 0o> is relevant, since the medium is horizontally homogeneous. It 
should be noted that with our definitions fi and fig are non-negative and we have 
< - 00 ^ 2;r. In this paper we are mainly interested in the question what 
happens with the reflection properties of the atmosphere if yu and yUo both approach 
zero, starting from values larger than zero. 

The intensity of the reflected radiation can be written in the form 

l'(fi,J^O, - 0o) = yUO^Ou.yUO, (/> - 0o)^O- (1) 

Here and hereafter the superscript t is used to indicate the top of the atmosphere and 
R(li,fio, 4> - (po) is the reflection function. The intensity at the top and the reflection 
function are both nonnegative. Once the reflection function has been obtained one 
readily finds the intensity of the reflected radiation for a parallel beam of incident 
radiation by using Eq. [T] while for multi-directional incident light an integration 
over all incident directions must be performed. 

3 Approaching the origin in the (/io,/i)-plane 

A function of one real variable may have two different limits at a point of the real 
number axis, namely a right-hand limit and a left-hand limit. For a function of 
two or more real variables there are many more possibilities to approach a point 
in the relevant multi-dimensional space and this may or ma y not correspon d to a 
number of different values to which the function approaches jCouran IT962h . This 
will be considered in the following sections for, respectively, the intensity of the 
reflected radiation, the reflection function and orders of scattering of both. 

3.1 Intensity of the reflected radiation 

Let us keep the azimuth difference fixed so that the functions /'(/j,yUo, if> - 0o) and 
RQj., fio, (f> — (po) are functions of the two variables, yu and /jq- We can approach 
the point yu = yUo = in various ways. This can be visualized by saying that we 
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can approach the origin, O, in the first quadrant of a Cartesian (jUo.A')"Coordinate 
system by following different paths (curves, including straight lines). Suppose we 
represent such a curve by means of a continuous function fi = g(fio) through O 
with a definite tangent at O (see Fig.lTJ. If we now approach the origin along this 
curve, the ratio g(/Jo)//Jo will tend to the slope of the tangent at O, which we denote 
as c. For a non-perpendicular tangent at O the value of c is finite and equals the 
right-hand derivative of g(/Jo) at O. Evidently, /jq and g(/Jo) are both non-negative, 
so that c is non-negative. If the curve is a straight line we simply have g(/Jo) = c/io- 
If we first let fi and then /jq tend to zero we have c = 0. If we first let fig and then /i 
approach zero we can treat this case separately or by letting c tend to infinity. 

As an illustration we may consider the special case of isotropic scattering in a 
semi-infinite homogeneous atmosphere with an albedo of s ingle scattering a. The 
phase function is then identically equal to one and we have jChandraseldQl 19501) 

/i(//,A<o,<^-0o) = — j — , (2) 

where the subscript 1 refers to the first order of scattering. Writing lim^.^|,^o for 
the limit if we approach O along a curve represented by /i = g(/io) we readily find 

Urn l\(^,fio,<p-<Po) = ., (3) 
tt-fo^o 4(c + i) 

which shows that the result depends on c, i.e. on the path that is taken to approach 
the origin. Consequently, there is a discontinuity for the intensity of the reflected 
radiation when the directions of incidence and reflection both beco me horizontal. 
This was called a peculiar discontinuity bv lHovenier & Stamll20o3) . since it looks 
at first glance rather surprising. 

In general, the scattering may be anisotropic, the atmosphere may be verti- 
cally inhomogeneous and its optical thickness may be finite with a reflecting or 
totally absorbing surface underneath the atmosphere. The intensity of the reflected 
ra diation in this general ca se for near-horizontal directions was also considered 
bv lHovenier & Stanil ( llOod ). They found a peculiar discontinuity in this intensity 
when ji and yUo both approach the origin in the (yUo,yu)-plane, which can be written 
as 

lim /V,y"o,0-<^o) = lim /',Ou,yUo,0 - ^o) 

Z'(cos(,^ - 0o))Fo, (4) 



4(c+ 1) 

where d is the albedo of single scattering at the top of the atmosphere and Z'(cos 0) 
is the phase function at the top of the atmosphere with scattering angle 0. Natu- 
rally, in Eq.|4]the same path must be followed for the two limits (see Fig. IB- Con- 
sequently, the following conclusions can be drawn for the intensity of the reflected 
radiation in the limit of /j and /jq both being zero: a) the optical thickness of the 
atmosphere and the reflection properties of the underlying surface are irrelevant, b) 
the values of the albedo of single scattering and the phase function need only to be 
known at the top of the atmosphere, c) orders of scattering higher than the first do 
not contribute, d) for any path with c unequal to infinity the azimuth dependence 
is proportional to the scattering angle dependence of the phase function at the top. 
The constant of proportionality becomes zero if c equals infinity. A similar state- 
ment was made bv lMinnaerll ( I1935D for a semi-infinite homogeneous atmosphere. 
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but he did not provide a coiTect proof and did not mention that the constant of 
proportionality depends on the way yu and /jq tend to zero. For more details about 
the behavior of the reflected intensity for dire ctions of incidence and re flection that 
both tend to horizontal directions we refer to lHovenier & Stanj ( l2006h . 



3.2 The reflection function 

Let us now consider what happens with the reflection function on approaching the 
origin in the (/Jci,/j)-plane. Combining Eqs.[T]and|4]gives 

lim ^o^(//,yuo,0-0o) = 

a' , 

hm fioRiipi,fio,(p-(po) = — — Z (cos(0 - (/>o)). (5) 

fi-m^o 4(c + 1) 

Using the principle of reciprocity jvan deHulstll980l) we have for all orders of 
scattering and their sum 

R(f^,fio, (p-(po) = R(f^o,f^, (p - (po)- (6) 

So, interchanging fj. and yuo in Eq.|5]and taking limits following the same paths as 
before we find 

lim iiR(jj,iJo,(t>-<l>Q) = 

hm iiRi(ii,iio,4)-(p(,) = — -Z (cos(0 - 0o)), (7) 

li-m^o 4(c + 1) 

since c in Eq. [5]had to be replaced by 1/c, i.e. the slope of g(jUo) at O with the 
positive ;u-axis. Since c is nonnegative, the limits in Eqs.|5]and|7]are bounded. So 
it follows from either one that 

lim nnoR{n,no, (j) - (po) = 0. (8) 

Since the right-hand side of this equation does not depend on c, this limit is the 
same for all curves represented by /j = g(jUo). More limits of this type, i.e. not 
depending on c, will be encountered further down in this paper. 
By adding Eqs.|5]and|7] we obtain 

lim (ji + iJ.o)R(ix,iiQ, (p-cpo) = 

d 

hm 0" +;Uo)/?iOu,yUo,0 - ^o) = — Z (cos(0 - 0o)). (9) 

fi.m^o 4 

This result has been reported for the spec ial case of a se mi-infinite, homogeneous 
atmosphere by seve ral authors (see e.g. ISobolevlfig?^ : [Mishchenko et ahi200^ : 
lKokhanovskvll200lh . but without a rigorous proof. 

Since the far right-hand side of Eq.|9]is positive we must have 



and 



lim i?i(;i,//o,<^ - 0o) = (10) 



lim i?Ou,yUo,0-0o) = oo, (11) 
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because zero or any finite number for these limits would be in conflict with Eq.[9] 
Consequently, we have proved that Ri(fi, fig, (^-0o) as well as R(ji,fio, (p-ipo) have a 
discontinuity if fi and /xq are both zero. The nature of these discontinuities is, how- 
ever, quite different from the peculiar discontinuities for the intensities discussed 
in Sect.im 

3.3 Orders of scattering 

The reflected intensity and the reflection function can be written as a sum (series) 
of nonnegative terms representing orders of scattering ( IvandeHulstI 19801) . In view 
of Eq.|4]we have for the n-lh order of scattering 

lim /;0L/,/io,'^-<^o) = O for n>l. (12) 

However, we cannot infer from Eqs. llOllTl what will happen with R„(pi,iio, <p - ipo) 
for « > 1 if // and fio both approach zero. The result might be zero, a finite positive 
number or infinity. However, some interesting properties for the higher orders of 
scattering of the reflection function can be obtained as follows. Equation [9] shows 
that 

lim (ji+fio)R„(jj,^Q,(l>-i/>o) = for n>l. (13) 
Furthermore, writing 

R{fi, fip, (p - ^o) _ (jU+yUo)/;Ou,jUo,(^-0o) ^^^^ 

Rl(fl,flO,<P - >Po) 0" +yt/o)i^lO",yUO,0 - ^o) 

and using Eq.|9]we find 

lim i?(p,/Jo,0-0o)//^i(p,A/o,0-0o) = 1, (15) 

since the right hand side of Eq.[9]is positive. Similarly, we obtain 

lim R„(ji,iuo,(l>-(f)Q)/Ri(ji,iuo,(l>-(f)Q) = 0, for n>l (16) 

The ratio of the reflection function to its first order term is an important function, 
since it is the correction factor to be applied to the easily computed first order term 
to obtain the reflection functi on. Numerica l studi es of this correction factor for 
homogeneous atmospheres by [van de HulstI Il980l) have shown that for isotropic 
scattering the correction factor has a maximum of 8.455, but that it becomes much 
smaller for nearly horizontal directions of incidence and reflection. He also found 
a similar behavior for anisotropic phase functions and mentioned Eg. 1151 

4 Examples 

To illustrate and check the results of the preceding section we will now discuss 
some simple examples. First for homogeneous and then for inhomogeneous atmo- 
spheres. 
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4.1 Homogeneous atmospheres 

Let us consider again tiie special case of isotropic scattering in a semi-infinite, 
homogeneous atmosphere with albedo of single scattering a. This implies that 
there is no azimuth dependence for the reflected radiation, so that we can omit 
(J - 00 in equations. The first two orders of scattering of the reflection function can 
readily be computed by analytic integrati on over o ptical depth (Hovenier 1 97 11) 
or by iteration of an invariance relation dAmbarzumiarji 19431 : iMishchenko et alj 
I2OOQ) . The results are as follows 



and 



Ri(p,^io) = ^ - (17) 



R2(u,iuo) = -(kUi) + ki^o)), (18) 



where 

k(M)=f,ln(l + l/fi). (19) 
For the sum over all orders we have dChandraseldi^l 19501) 

R(M,tio) = ,, " H(ji)H(Mol (20) 

where H(ji) is a well-known function, depending on a, with H{0) = 1 ( lBusbridg3 
1 1960 ). Equations 1 171201 show that the reflection function and its first two orders 
of scattering tend to infinity if fi and fjg both approach zero. Consequently, this 
also holds for the multiple scattering component of the reflection function. Clearly 
Has. 1 171201 are in agreement with Eqs. 1511 Il and ll3l76l since k(0) = 0. Figure |2] 
shows Ri(ji,iU(,) and R(ji,fio) as functions of yu if a = 1 and fi = yUo- Both functions 
are seen to be strongly increasing when the directions of incidence and reflection 
tend to horizontal directions. 

It follows from Eqs.[T7]and[20]that the correction factor 

R(ji,fio)/R,(fi,fio) = mM)H(Mo). (21) 

This factor varies between 1 (if yu = /jq =0 and a is arbitrary) and 8.455 (if 
ji = jiQ = I and a = 1). Fig.[2]also shows this correction factor in case /j = fig. 

A lesson to be learned from this simple case is that one should not assume that 
the discontinuity of the reflection function for horizontal directions of incidence 
and reflection will always disappear upon integration. Indeed we find from Eq.|17| 



Jo 



dfioR,(fi,fio) = ^Ml + i-/l^), (22) 



which tends to infinity if fi tends to 0. So this approach to infinity must also hold 
if the reflection function is integrated in the same way, since the sum of all orders 
of the reflection function cannot be smaller than the first order only. Another way 
to prove this can be obtained from the definition of the //-function written in the 
form ^ 

r d^loR(^l,^lo) = ^(H(n)-i) (23) 

Jo 2yU 

and the behavior of the //-function when fi approaches zero jvan deHulstll 19801) . 
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Explicit expressions for the reflection function of a semi-infinite or finite ho- 
mogeneous atmosphere on top of a blac lc surface in terms of functions H(ji), X{ju) 
and Y(p) have been pubhshed (see e.g. IChandrasekhaill 19501) for phase functions 
that can be written as a sum of a few Legendre polynomials, including Rayleigh 
scattering. Since H(0) = 1, X{0) = 1 and Y(0) = (Busbridge 1960) it can readily 
be verified that in all these cases Eos.lSHTTIand Eq.[T5]are valid. 



4.2 Inhomogeneous atmospheres 

In numerical calculations a vertically inhomogeneous atmosphere is often modeled 
as a stack of homogeneous layers. Then only the albedo of single scattering and 
the phase function of the top layer are relevant for the reflected intensity when /i 
and fio both approach zero. When these values are substituted in Eqs.|4l[5]|7l and[9l 
Eqs. 141161 are also valid for a stack of homogeneous layers. 

Another example of scattering in an inhomogeneous atmosphere is provided by 
isotropic scattering in a semi-infinite atmosphere with albedo of single scattering 
given by 

a(T) = exp(-T), (24) 

where r is the optical depth me asured from the top of th e atmosphere downwards. 
In this case we have for /jq > faovenier & Stamll2006r) 

l\(p,fio) = ^ . ^ ^° (25) 

so that 

R>{M,f^o) = \ . ^ \ -. (26) 

It should be noted that the denominators in the last two equation s are not simply the 
sum of /J and /uq. Eqs. |25tj26l and the expressions reported bv lYanovitskiil ( ll997D 
for the total intensity of the reflected radiation and the reflection function are in 
agreement with Eqs. 141161 So the latter equations provide a useful check. 



5 Related functions 

In addition to the reflection function several alternative functions for describing 
the reflective properties of a plane-parallel atmosphere are found in the literature. 
If the diflFerence with the reflection function is only a constant factor all preceding 
equations must be translated by simply taking this constant into account. But the 
difference may also involve a function of /j and/or /jq and then the translation is 
less trivial. In such cases one can derive expressions for the alternative function by 
starting with an expression for the intensity of the reflected light in terms of the al- 
ternative function and then working along similar lines as we did for the reflection 
function or directly use relations for the reflection function to get corresponding 
res ults for the alternative function. We give some examples. 

IChandrasekhajh950l) defined what he called the scattering function by writing 

Fo 

I'(ji,Ho,^ - (po) = —S(ji,4i;fio,<po), (27) 
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so that in view of Eq.[T]the scattering function is related to the reflection function 
as follows 

S(ji,(p;fio,(f>o) = 4fifioR{ii,iio,(l>-(/>o). (28) 
It follows immediately from this equation and Eq.[8]that 

lim S(p,<p-fio,<l>o) = 0. (29) 

So the same must be true for all orders of scattering, since we are dealing here 
with a sum of nonnegative numbers. Consequently, the orders of scattering of the 
scattering function, nor the function itself, tend to infinity if we approach the origin 
of the (//o,y")-plarie along a path given by yu = g(j-io)- Using Eq.|9]we find 

lim (l/yu+ \/fio)S(fi,'l>;fJo,(f>o) = 

lim (1/yU + l//io)5iO",(^;A'o,<^o) = a'Z'{cos{^ - (po)), (30) 

which again shows the special role played by single scattering. 

Some authors still use the scattering function, but the reflection function is 
more w idely used today. Each one has certain advantages and disadvantages. 

defined several alternatives for the reflection function, most of which 
difl^er only by a constant factor from the reflection function. But his so-called bidi- 
rectional reflection function r(ji,fio, (p - (f>o) has the same peculiar discontinuity as 
our f(ji,fio,(/>- (pojIinFo). 



6 Azimuth dependence 

So far we have kept the azimuth difl^erence (p - <po fixed and considered what hap- 
pens when not only the incident radiation tends to the horizontal direction in the 
(^o-plane but also the reflected radiation tends to the horizontal direction in the (bi- 
plane. More generally, the azimuthal angles may change when ^ and tend to 
zero, i.e. {pi, (/>) may tend to (0, 0) and (pio, (pg) to (0, (po) in some way. Clearly, the 
discontinuities for a plane-parallel atmosphere will then remain and we must only 
replace cos(0 - <Pq) by cos(0 - <po) in Eqs.|4l|5]|7]and[9] Such situations generally 
occur in spectrophotometry of planets and satellites when one moves the line of 
sight to the intensity poles (jj = fio = 0), along the limb (pi = 0), the terminator 
(/io = 0) or an intermediate intensity meridian (p = CyUo) with finite c > 0. This 
must be realized in c omputat ions when the outer layers are modeled as a plane- 
parallel atmosphere teovenier & Stam 2006). 

The azimuth dependence of intensities and reflection functions are often han- 
dled by making Fourier series expansions, e.g. by writing for an arbitrary function 

f(M,^io,'p- M = /V.A'o) + 2E,:^,/"V,A'o)cos(m(0 - 0o)), (31) 

where the upper index m denotes the Fourier index. In this way we readily find 
expressions for each Fourier component separately from the equations given in 
preceding sections. For instance, writing Z'(0, 0,ip - (po) for the phase function at 
the top of the atmosphere in case yU = /io = 0, we obtain for in > using Eqs. |5] 
and|7] 

lim //o R"'(p,f^o) = — ^Z''"(0, 0) (32) 
ti-m^o 4(c + 1) 
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and 



lim u i?"V,y"o) = -^^Z'"'(0, 0). (33) 
f'-m^o 4(c + 1) 



Adding the last two equations gives 

lim 0" + t^o)R"'(p,l^o) = tZ'-'"(0, 0). (34) 

The azimuthal average of the phase function is positive, since we have assumed 
that the phase function itself is positive. Hence, the m = component of the phase 
function is positive and it follows from Eq.|34|that 

lim R°(n,fio) = oo. (35) 

However, Z''"(0, 0) can be positive, negative or zero if m > 0. Consequently, if fi 
as well as yUo tend to zero the limit of R"'(ji,iU(,) may be plus infinity, minus infinity 
or it may be impossible to determine this limit from Eq.[34l This may be illustrated 
by considering the following simple cases. 

(i) The phase function 

Z'(cos 0) = 1 - - cos (36) 

has 

Z'-°(0, 0) = 1 (37) 

and 

Z'i(0,0) = -1/4. (38) 

(ii) Rayleigh's phase function, i.e. 

Z'(cos0) = ^(1 +cos2 0) (39) 

has 

Z'''(0,0) = 9/8 (40) 

Z''(0,0) = (41) 

Z''^(0,0) = 3/16. (42) 

Using the explicit expressions l lChandrasekhai|[l95Cl) for the reflected intensity of 
a homogeneous semi-infinite atmosphere with a = 1 for these two phase functions 
one can easily find that Eqs. |32l[34l are satisfied. 



7 Approximations 

Accurate computations of the reflection function of a plane-parallel atmosphere 
are generally not easy and also laborious, especially when results are needed for 
many model parameters. Therefore, it is not surprising that a variety of approxima- 
tive f ormulae have been proposed ( Sobolev 1975; van de Hulst 1980; Kokhanovskvl 
I2OOI ). Here we wish to point out that at least some of these are not realistic for 
directions of incidence and reflection that are both nearly horizontal. 
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For instance, according to the popular Lambert's reflection law the reflection 
function is a constant not larger than one for all directions of incidence and reflec- 
tion. So this reflection function does not obey Eq.[9]and can never tend to infinity, 
which is in conflict with Eg . 1 1 1 1 

Another example is the so-called "rapid-guess formula" of lvan deHulstlfT980l) 
for a non-absorbing homogeneous semi-infinite atmosphere 

^V^A-o) = 1 + P(l - 3a</2)(1 - 3/io/2), (43) 

where the upper index, 0, refers to the azimuth independent term and p is a finite 
constant. This formula clearly violates Eq.l35land. therefore, may give unaccept- 
able errors for nearly horizontal directions of incidence and reflection. To show this 
we consider isotropic scattering for which the value p = 0.4 was recommended by 
van de Hulst (1980). Using this value in the rapid-guess formula with ^ = /io = 0. 1 
gives for the reflection function 1.289, which is 33.72 % too low, since accurate 
multiple scattering calculations give 1.94485. For smaller values of p and po the 
errors in the reflection function are still larger. Integration of the reflection func- 
tion as in Eg. 123 1 gives for /j = 0.1 the accurate value of 1.23675, but Eg. 143 1 gives 
1.0850 which is still 12.27 % too low. 

Finally, we mention a formula that is often used for a homogeneous atmosphere 
above a black surface when the optical thickness b is z small positive number, 
namely 

ab 

R(M,I^o, <P-<Po)= Z(cos 0), (44) 

4ppo 

where 

cos = -ppo + -p^)(l -pl)cos{if> - 0o)- (45) 

Although Eq.|44]does not violate Eq.[TT]the approach of this approximate reflection 
function to infinity if p and po both tend to zero is apparently not correct, since 
Eas.l5]and l7l9l are not satisfied. 

Consequently, if p and po are very small one should be careful with using 
approximations for the reflection function. Fortunately, the results of the preceding 
sections suggest th at it may then be sufficient to compute only a few orders of 
scattering (see e.g. lHovenieil[T97ll: Ivan de Hulstlll98of : iMishchenko et alj|2006l) 
instead of a more laborious complete multiple scattering calculation. This holds 
in particular when the optical thickness and/or the albedo of single scattering is 
not large. For example, for isotropic scattering in a homogeneous semi-infinite 
atmosphere with a = 0.4 and p = po = OA the sum of the first two orders of 
scattering of the reflection function is only 1.6 % too low and the sum of the first 
three orders even less than 0.34 %. 

8 Concluding remarks 

Numerous complicated equations occur in the theory of multiple light scattering in 
hom ogeneous and inhomogeneous plane-parallel atmospheres (see e.g. ChandrasekhaJ 
1950 - Sobolev 19 75; van de Hulst 1980; Yanovitskij 1997; Hovenier et al. 2004; 
Mishchenko et aljfeooeh . Since there is always a possibility that printed equations 
contain errors and their derivations are not always given, it is useful to have simple 
checks available like letting p and po approach zero and comparing the results with 
expressions in this paper. 
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On performing model computations it is usually very helpful to know and un- 
dersta nd what happ ens with the reflection function or a related function in limiting 
cases ( llrvinein"983h . This kind of knowledge is provided in this paper even for a 
complicated model of an inhomogeneous atmosphere with an arbitrarily reflect- 
ing surface underneath. In particular one should be prudent in the proximity of 
discontinuities like those presented in this paper. 

We have shown that a discontinuity exists for the reflection function which 
may hamper interpolation and extrapolation. This problem may be by-passed by 
multiplying the reflection function by a simple function of fi a nd fip, like n + no (c f. 
Eq. |9]l, before performing the interpolation or extrapolation jKnibbe et al.ll200oh . 
This is illustrated in Fig.[2]by the curve marked R/R] since in the case considered 
this curve equals 8yui?(jU,/j). 

We have also shown that great care should be exercised with using approxi- 
mative formulae for the reflection function, since they may lead to large errors for 
nearly horizontal directions. This holds, for instance, for the "rapid guess formula" 
and the very popular Lambert re flection law, which is often used, e.g. for cloudy 
atmospheres of planets (see e.g. iKokhanovskvl 1200 6'). Although approximative 
formulae exist that do not give large errors for nearly horizontal directions, more 
accurate results are possibly obtained by computing a few orders of scattering for 
such directions. 
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Figure 1: A point P approaches the origin O of a Cartesian coordinate system along a 
curve (solid line) represented by p = gijio), which has a slope c at O with respect to the 
positive juo-axis. The tangent of the curve at O has also been drawn (dashed line). 
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Figure 2: The reflection function (R) of a non-absorbing, homogeneous, semi-infinite 
atmosphere with isotropic scattering in case jj = yuo is plotted as a function of 
/i. Also shown are the contribution due to first order scattering (Ri) and the ratio 
R(ji,fj)/Ri{fi,lj), which in this case equals 8yu/?(ju,//). 
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